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A study of the breakdown of the quasi-static approximation at high densities and its 
effect on the helium-like Ka complex of nickel, iron, and calcium 

Justin OelgoetzQ Christopher J. Fontes, and Hong Lin Zhang 
Applied Physics Division, Los Alamos National Laboratory, Los Alamos, NM 87545 

Anil K. Pradhan 

Astronomy Department, The Ohio State University, 140 W. 18th Avenue, Columbus, OH 43210 

Recent work to include R-Matrix data within a larger model comprised mostly of distorted-wave 
and plane-wave-Born data has resulted in the General Spectral Modeling (GSM) code. It employs 
a quasi-static approximation, a standard, low-density methodology that assumes the ionization 
balance is separable from a determination of the excited-state populations that give rise to the 
spectra. GSM further allows for some states to be treated statistically as contributions to effective 
rates, instead of being included explicitly in the kinetics model. While these two approximations 
are known to be valid at low densities, this work investigates using such methods to model high- 
density, non-LTE emission spectra and determines at what point the approximations break down by 
comparing to spectra produced by the Los Alamos National Laboratory code ATOMIC which makes 
no such approximations. As both approximations are used by other astrophysical and low-density 
modeling codes, the results should be of broad interest. He-like Ka emission spectra are presented 
for three elements, Ni, Fe, and Ca, in order to gauge the effect of both the statistical methods and 
the ground- state- only, quasi-static approximation employed in GSM. This work confirms that at 
and above the temperature of maximum abundance of the He-like ionization stage, the range of 
validity for both approximations is sufficient for modeling the low- and moderate-density regimes 
one typically finds in astrophysical and magnetically confined fusion plasmas. However, a breakdown 
does occur for sufficiently high densities; we obtain quantitative limits that are significantly higher 
than previous works. Additionally, this work demonstrates that, while the range of validity for both 
approximations is sufficient to accurately predict the density- dependent quenching of the z line, the 
approximations begin to break down at higher densities. Thus these approximations should be used 
with greater care when modeling high-density plasmas such as those found in laser-driven inertial 
confinement fusion and electromagnetic pinch devices. 



I. INTRODUCTION 



The subject of Ka emission lines arising from helium- 
like ions has been a topic in the literature for the bet- 
ter part of a century The work that the present ef- 
fort builds on begins with Gabriel and Jordan's work on 
the Sun 0, H| • There are four Ka lines arising from the 
helium-like ionization species, each involving decay from 
an excited state to the ground state. The 'w' line repre- 
sents a dipole allowed transition arising from the decay 
of the \s2p 1 P% state. Because of mixing, the \s2p 3 P% 
state can also decay via a dipole allowed transition, giving 
rise to the intercombination line 'y'. The ls2s 3 Si state 
can decay via a relativistic magnetic-dipole transition, 
giving rise to the forbidden line, 'z'. Lastly an additional 
forbidden line, 'x', is formed from the \s2p 3 Pg state de- 
caying via a magnetic-quadrupolc transition. Due to the 
sometimes small energy separation between the \s2p 3 P° 
and \s2p 3 P% states, the x and y lines can not be inde- 
pendently resolved for all elements or plasma conditions. 
The ls2s 1 S'o state must also be included explicitly in 
the determination of the excited-state populations as it 
can decay via a two-photon transition which has an im- 



portant effect on the model, as does the presence of the 
ls2p 3 Pg state @,||. 

From these four lines, line ratios have been developed 
in order to diagnose certain conditions about the corre- 
sponding plasma. One of those line ratios, R, is defined 
as 



R = 



J(z) 



f(x)+/(y) 



(i) 
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The intensities, /, of the non-dipole allowed lines (x and 
z) and the intercombination line (y) vary in a similar 
manner as a function of temperature. However, electron 
collisions can transfer population from the ls2s 3 Si state 
to the \s2p 3 Pq 1 2 states, thereby quenching the z line 
and increasing the intensity of the x and y lines. There- 
fore the ratio R can be used as a diagnostic of electron 
density. 

Building on this earlier work, the study of Ka emission 
from He- like ions was extended to heavier elements by in- 
cluding the effect of diclcctronic satellite lines @ , and to 
a broader range of elements and transient conditions in 
the seminal works of Mewe and Schrijver 0, [|| . Improved 
data and improved models continued to be brought to 
bear as they became available 0, H, H, floL ll[ . Photoion- 
ized media have also been investi gate d [12 1 . Lastly, more 
recent Breit-Pauli R-Matrix [l3l Il4j data were used to 

refine the quantities of interest even further mill. Ad- 
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ditionally the Ka lines are also observed in laboratory 
plasmas [J, EE GJ, EI US HH, H3, HI, and these same 
methods have been successfully used to model and inter- 
pret the spectra of magnetically confined fusion plasmas 
|19l . I20I , [2ll . |22| . While the above narrative and cita- 
tions are only a small portion of the work done in this 
field since these lines were initially observed, they are 
the platform upon which the General Spectral Modeling 
(GSM) code [2J], one of the codes used in the present 
work, is based. 

It should be noted that much of the previous work ex- 
plicitly assumed that, for the purposes of finding the pop- 
ulations of the He-like excited states, the populations of 
the adjacent ionization stages are entirely in the ground 
state. While this particular type of quasi-static approx- 
imation (the ground-state-only, quasi-static approxima- 
tion) is clearly valid at low densities, it breaks down at 
higher densities as electron-impact ionization and recom- 
bination out of excited states become important mecha- 
nisms for populating the excited states. One could en- 
vision extending this treatment by including not only 
the ground state in the set of states with which the ex- 
cited states are in equilibrium, but also select metastable 
states. The effect of expanding the quasi-static treat- 
ment employed in this way is not explored in the present 
work, as the vast majority of papers on astrophysical 
He-like Ka emission make these same approximations. 
Additionally many of the above works treat many of the 
excited states as cascade corrections (see, for example, 
@) to direct rate coefficients. The effect of such a statis- 
tical approximation is explored in the present work. 

From a more general perspective, expanding beyond 
the specific case of He- like ions described above, the trend 
over time has been to incorporate more accurate atomic 
data in plasma kinetics and spectral modeling. Often, 
these better atomic data are difficult or impossible to cal- 
culate for all processes between all levels. Furthermore, 
if the atomic data are calculated using a method that 
accounts for resonance structure, such as the R-Matrix 
method, it is much more time consuming to calculate 
rate coefficients from the more complex cross sections 
than from smooth cross sections. Reducing the compu- 
tational requirements makes the quasi-static and statis- 
tical approximations appealing. A separate, but related, 
consideration is that as the plasma transitions from non- 
local-thcrmodynamic-equilibrium (non-LTE) conditions 
to the high density regime, an accurate modeling ap- 
proach should produce the appropriate LTE limit. How- 
ever, each of the two approximations mentioned above 
typically preclude a transition to the correct LTE pop- 
ulation distribution in the high-density limit, regardless 
of the quality of the atomic data that is used. The flexi- 
bility of GSM allows for an investigation of the interplay 
between these various considerations when attempting to 
model plasmas over a broad range of conditions. 

The primary interest of this work is benchmarking and 
examining the range of validity for the ground-state-only, 
quasi-static approximation as well as the statistical treat- 



ment of excited states in the context of He-like ions, as 
these approximations are employed by most of the works 
referenced above. Many of the connections between ex- 
cited states are neglected by these two approximations; 
thus the constraints on the system are not sufficient to 
produce the correct high-density, local thermodynamic 
equilibrium (LTE) limit. Using these approximations re- 
sults in a significant speed up however, and allows one 
to include the effects of many more states with the same 
computing power. Thus the applicability of these approx- 
imations to high-density, non-LTE plasmas such as those 
found in vacuum sparks [l7| , electromagnetic pinch type 
devices [23l . |25| and laser produced, inertial confinement 
fusion plasmas [HJ is of great interest. To gauge the effect 
of these approximations, the spectra produced by GSM 
are compared to the results from another code, ATOMIC 
[2^ |. ATOMIC docs not employ cither the quasi-static 
approximation or statistical methods for treating excited 
states; instead every fine-structure level in the model is 
treated identically, such that all included processes be- 
tween all levels are considered explicitly It should be 
noted that the present work is not the first compari- 
son of a code based on the quasi-static approximation to 
ATOMIC; comparisons between the code ADAS [27| and 
ATOMIC have been previously considered [2q . However, 
the primary aim of the present work is to investigate the 
breakdown of the approximations involved and to provide 
quantitative limits for the breakdown, while the earlier 
work considered the effects of differing atomic data sets 
on plasma kinetics modeling. 

The secondary aim of this work is to confirm that using 
a statistical treatment and the ground-state-only, quasi- 
static approximation are adequate for modeling emission 
spectra in the low- and moderate-density regimes encoun- 
tered in magnetic confinement fusion devices and most 
astrophysical plasmas. As the line ratio R is useful for de- 
termining the density of the emitting plasma, and many 
previous calculations of R employ both a statistical treat- 
ment of excited states and the ground-state-only, quasi- 
static approximation, we also confirm that the methods 
used for such calculations arc indeed adequate. The test 
cases presented in this work arc for Ni, Fe, and Ca as 
they arc directly relevant to astrophysical and laboratory 
plasmas. 



II. THEORY 

Two distinct atomic kinetics codes, ATOMIC (2(| and 
GSM [HI, have been used to explore the same plasmas 
in an effort to determine the range of validity of the 
ground-state-only, quasi-static approximation and statis- 
tical treatment of excited states for the three elements 
considered in this work. ATOMIC has evolved within the 
context of modeling high-density plasmas where assump- 
tions like the quasi-static approximation can be prob- 
lematic. As previously mentioned, GSM employs the 
ground-state-only, quasi-static approximation, and has 
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the ability to treat states statistically. While GSM en- 
ables the use of more accurate atomic parameters, such 
as R-matrix data, the methods used by GSM are an ap- 
proximation to the more general approach considered in 
ATOMIC; thus it is convenient to start with a descrip- 
tion of the theory employed by the ATOMIC code and 
then move on to describing GSM. 



s(hu) =^1(1,3 - hk)hv^± r-V-r' ) (4) 

where nij and Tj are the ion mass and temperature re- 
spectively, and AEjk is the energy of the line produced 
by radiative decay from level j to k. 



ATOMIC 



B. GSM 



ATOMIC [26| is a general kinetics code which solves 
the collisional-radiative atomic equations given by 



dN, 



id 



dt 



Ne £ {Nuq^e) - N^q^e)) 
£Ni + x,iCi +hi ^i d {e) - Ni, j £Ci, j ^i- 1 , i (i) 

i i 



Ni 



-l,i a l-l.i^Lj( £ ) 



l+l,i 



(i<j) 



A, 



Ni 



£ Ni- lt iPi- lti ->ij (g) 



i I i(i>j) 



i(i>j) 

£ Ni+^iRf+ij^j - Nij £ R 



AI 



(2) 



One can see from Eq. © that all states are treated ex- 
plicitly whether they be bound or autoionizing. Rate 
coefficients for electron- impact excitation (q), collisional 
ionization (C), dielectronic capture (D BC ), radiative re- 
combination (a RR ), three-body recombination (/3), ra- 
diative decay (A), and autoionization (R Al ) are included 
between all possible fine-structure levels in the model. 
A variable describing the electron energy distribution is 
denoted by e; as all results presented in this work are 
for thermal systems, e can be taken to be the electron 
temperature. As radiation fields are not considered in 
this work, processes depending on them have been omit- 
ted from the discussion but are in general included by 
ATOMIC. This method is guaranteed to go to the cor- 
rect LTE limit as the electron density increases since the 
rate coefficients for the inverse processes are calculated 
assuming detailed-balance with those of the forward pro- 
cesses. 

Once the populations (Nij) have been calculated for 
each level (j) of every ionization stage (I), emission spec- 
tra (S) are determined by finding the intensity (I) of each 
line and then applying a Doppler-broadened line profile 
to it using the equations 



GSM [24| is based on the quasi-static approximation 
(see, for example, [28j |) as implemented by much of the 
previous work on modeling He-like spectra [1, 0, H, HE 
ll6f| . As such, GSM assumes that the excited states of 
a particular ionization stage are always in instantaneous 
equilibrium with the adjacent ionization stages, which 
are assumed to be entirely in the ground state. The as- 
sumption is valid at low densities since the excitation and 
radiative decay rates are generally higher than ioniza- 
tion and recombination rates that determine ionization 
balance. 

The result of employing the ground-state-only, quasi- 
static approximation is that the first step in any calcu- 
lation is to obtain the ionization fractions by solving the 
following set of coupled equations 



dXi 
dt 



1(1, j -+l,k) = N Lj A j ^ k , 



(3) 



N e (Xi +1 at +1 .+i(e) + AViG-i^O?)) 
+Nl(X l+1 p l+1 ^ l (s) - XiA^-i(e)) 
-AW e («,->/_!(£) + CU I+ i(e)) , (5) 

where Xi is the total population in the I th ionization 
stage, N e is the electron number density, C is a bulk col- 
lisional ionization rate coefficient, (3 is a bulk three-body 
recombination rate coefficient, and a is a bulk recom- 
bination rate coefficient (which includes radiative and 
dielectronic recombination). Again, terms involving a 
radiation field, such as photoionization or stimulated re- 
combination, have been omitted as this work assumed a 
collision-dominated plasma for which one can neglect the 
radiation field. The bulk rate coefficients arc typically 
calculated by summing over the individual pathways in- 
volved, or from literature sources that have performed 
such operations (e.g. [2{|), and the Xi are obtained by 
solving the resulting set of coupled equations. However, 
in this study the solutions to these equations are taken 
to be those found by summing over the level populations 
(Ni.j) of each ionization stage as determined from an ex- 
plicit ATOMIC calculation (sec Eq. © and Fig. [TJ). 
This choice removes the possibility of any discrepancy 
that might be caused by inconsistencies in the ionization- 
balance data employed by each code. 

GSM offers the additional option to treat a portion of 
the states in a statistical manner, as conduits or inter- 
mediate pathways that arc 1 involved in the calculation of 
the rate coefficients that are used in a reduced system 
of coupled equations. The populations of these conduit 
states are not computed, but all other levels are treated 
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explicitly. The result is a set of effective rate coefficients 
which, along with the total populations in each ioniza- 
tion stage, are then used to determine the excited-state 
populations by solving the system of coupled equations 
given by 



1,3 



(It 



where sums correspond not only to excitation followed 
by radiative cascade, but also correspond to capture into 
an autoionizing state followed by autoionization, or some 
combination of autoionization and radiative decay which 
produces the final explicit state k. Similar expressions 
hold for other effective rate coefficients. Once the effec- 
tive rate coefficients are used to calculate the level popu- 
lations via Eq. ([5]), these populations are used to produce 
spectra using Eqs. (JSJ) and @, just as in the ATOMIC 
calculations. 



X, 



lie) - N u J2<l- 



u+i(e) III. COMPUTATIONS 



E 

i(i>j) 



X, = 



i 

Z 



i(i<j) 



? AI-off 



,3 > 



where the terms are defined identically to the case for 
ATOMIC, except that the superscript 'eff ' denotes an ef- 
fective rate coefficient, which includes contributions from 
levels that do not appear explicitly. Just as in Eqs. © 
and iJSJl, Eq. © omits terms that involve the radiation 
field since it is neglected in this work. 

These effective rate coefficients are calculated by sum- 
ming over the direct and all indirect paths involving 
states being treated statistically. The indirect paths 
involve the use of the collisionless transition matrix 
(CTM), Tto-jj, which can be thought of as the proba- 
bility that an ion in statistical state m will end up in 
state j, via spontaneous processes. Thus radiative decay 
and autoionization are the only processes included in cal- 
culating the CTM. It is calculated from these rates using 
the recursive relation 

rp _ \ ■* ry_>fc 

*-i-*3 — A , dAI Lk ^3 



E^ + E^, 

{E i >E k >E j ) I m 

r 



E^ + E^, 



where Ti—>k is either a radiative decay rate or an autoion- 
ization rate. The CTM is used in the calculation of all the 
effective rate coefficients. For example, effective electron- 
impact excitation and de-excitation rate coefficients are 
calculated as 



(E i >E j ,E i >E k ,E i >0) 



In order to simplify the comparisons, GSM was run 
with the ionization-balance data calculated by ATOMIC, 
as mentioned earlier in the discussion following Eq. ((5]). 
While the ionization balance is a function of both electron 
temperature and density, the temperatures considered for 
each element were chosen relative to the low-density limit 

(6) of the ionization balance (see Fig. []}. The first temper- 
ature (5.17 x I0 3 eV (6.0 x I0 7 K) for Ni, 3.02 x 10 3 eV 

(7) (3.5 x I0 7 K) for Fe, and 1.29 x I0 3 eV (1.5 x 10 7 K) 
for Ca) is near the temperature of maximum abundance 
of the He-like ionization stage. The second tempera- 
ture (1.72 x 10 4 eV (2.0 x 10 s K) for Ni, 8.62 x 10 3 eV 
(1.0 x 10 8 K) for Fe, and 3.45 x 10 3 eV (4.0 x 10 7 K) 
for Ca) is near the temperature of maximum abundance 
for the H-like ionization stage. Lastly the third temper- 
ature (4.31 x 10 4 eV (5.0 x 10 s K) for Ni, 2.59 x 10 4 eV 
(3.0 x 10 8 K) for Fe, and 6.89 x 10 3 eV (8.0 x 10 7 K) 
for Ca) is significantly above the previous two tempera- 
tures and considers the case for which the bare nucleus is 
dominant. Furthermore, both ATOMIC and GSM were 
run with the same set of fundamental atomic data, calcu- 
lated with the Los Alamos suite of atomic physics codes 
[30l . l3~fl ] . The CATS code was used to calculate all wave 
functions, energy levels, dipole allowed radiative decay 
rates, and plane-wave-Born electron-impact excitation 
cross sections arising from bound and autoionizing levels 
of the configurations nl, Isnl, 2lnl', ls 2 nl, \s2lnl 1 ', and 
IsZlnl' with n < 10 and I < g. The output of CATS 
was then used by the GIPPER code to calculate au- 
toionization rates and photoionization cross sections in 
the distorted-wave approximation, as well as collisional 
ionization cross sections using a scaled hydrogenic ap- 
proximation which accurately reproduces distorted-wave 
results for the ionization stages of interest. The ACE 
code was used to obtain distorted-wave, electron-impact 
excitation cross sections for all transitions out of the low- 
est seven levels of the He-like ionization stage, as well as 
out of the three levels in the ls 2 2l complex of the Li-like 
stage. Lastly, non-dipole radiative decay rates for the 
three elements under consideration, i.e. the magnctic- 
dipolc, magnctic-quadrupole, and the two-photon rates 
mentioned previously, were obtained from Mewe and 

D^(e)Ti^ k , (9) Schrijver [|. 

As described above, two approximations are consid- 
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FIG. 1: (Color online) Ionization-balance data determined by ATOMIC for the three elements under consideration at an 
electron density of iV e = 10 10 cm -3 . The three temperatures considered for each of the three elements have been indicated. 
The first temperature (indicated by the arrow labeled Ti) was chosen to be near the temperature of maximum abundance for 
the He- like ionization stage (the solid red line), the second temperature (indicated by the arrow labeled T2) was chosen to be 
near the temperature of maximum abundance for the H-like ionization stage (the dashed green line) , and the third temperature 
(indicated by the arrow labeled T3) was chosen such that a significant fraction of the ionization balance consisted of the bare 
nucleus (the blue dash-dot-dot line). 



ered for the GSM calculations: (1) some states can be 
treated statistically, and (2) the excited states are in in- 
stantaneous equilibrium with just the ground states of 
the appropriate ionization stages. In order to test the 
validity of the statistical treatment of certain states, two 
GSM models were considered. The first model, called 
GSM7, treats the lowest seven levels of the He-like ion- 
ization stage, along with the ground state of the H-likc 
ionization stage, and the levels arising from the ls 2 2l 
and ls2lnl' configurations in the Li-like ionization stage 
as explicit. The second model, called GSMF, considers 
all levels to be explicit. Thus, this full model includes and 
accounts for collisions among the excited states within an 
ionization stage. It should be noted that due to the sec- 
ond approximation, this full model neglects population 
mechanisms into these excited states from levels other 
than the ground state in the adjacent ionization stages. 
It does, however, include recombination and ionization 
from these excited states into all the levels of the adja- 



cent ionization stages as possible mechanisms for popu- 
lation loss via the effective rate coefficients. Whenever 
results obtained from the GSMF model differ from those 
computed with GSM7 model, it is an indication that the 
excited states can no longer be treated statistically. 

Similarly, the validity of the ground-state-only, quasi- 
static approximation is tested by comparing results from 
the GSMF model to the output of the ATOMIC code 
(which is referred to as ATOMIC). Differences between 
the GSMF and ATOMIC models indicate that the sepa- 
ration of the ionization-balance and excited-state popu- 
lation calculations is not valid. 



IV. RESULTS AND DISCUSSION 

The first element to be considered is nickel (see Fig. [5]) . 
The top two rows of Fig. [5] show emissivities generated 
from the three models that arc in excellent agreement. 
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FIG. 2: (Color online) Ni emission spectra at T e =5.17x 10 3 , 1.72 x 10 4 , and 4.31 x 10 4 eV (6.0 x 10 7 , 2.0 x 10 8 , and 5.0 x 10 8 K) 
for the three models considered at four different electron densities (top to bottom): 10 10 , 10 19 , 10 21 , and 10 24 cm -3 . The solid 
line (green) was calculated with the ATOMIC model, dashed (blue) with the GSM7 model, and the dot-dot-dashed (red) with 
the GSMF model. 



The top row of spectra was calculated at a relatively 
low electron density (10 10 cm -3 ), typical of astrophys- 
ical sources. The second row of spectra was calculated 
at a much higher density (10 19 cm -3 ), but one at which 
all three models still agree. This agreement defines a 
range of validity which persists among all three temper- 
atures. The unlabeled features are dielectronic satellite 
lines. They are formed only at low temperatures, and 
become weaker at higher densities, most likely due to 
collisional quenching. By the third row (10 21 cm -3 ), the 
GSM7 model has begun to disagree with the ATOMIC 
and GSMF models, which is an indication that the sta- 



tistical treatment of excited states is no longer valid. In 
the bottom row (10 24 cm -3 ) all three models are vastly 
different at the lowest temperature, an indication that 
both the ground-state-only, quasi-static approximation 
and the statistical treatment arc no longer valid. How- 
ever, the ground-statc-only, quasi-static approximation 
remains valid at the highest temperature. This trend 
has a simple explanation: at higher temperatures, ion- 
ization out of the Li-like ionization stage plays less of 
a role in He-like spectral formation. Even though Li- 
like metastable states are important fractions of the Li- 
like ionization stage population, the effect of collisional 
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FIG. 3: (Color online) Fe emission spectra at T e =3.02 x 10 3 , 8.62 x 10 3 , and 2.59 x 10 4 eV (3.5 x 10 7 , 1.0 x 10 s , and 3.0 x 10 8 K) 
for the three models considered at four different electron densities (top to bottom): 10 10 , 10 19 , 10 21 , and 10 24 cm~ 3 . The solid 
line (green) was calculated with the ATOMIC model, dashed (blue) with the GSM7 model, and the dot-dot-dashed (red) with 
the GSMF model. 



ionization from them into the He-like ionization stage is 
not as large of a factor due to the overall decrease in 
the population of the Li-like ionization stage. In addi- 
tion, the populations of the He-like and H-likc ionization 
stages remain dominated by their ground state values. 
Thus the ground-state-only, quasi-static approximation 
remains valid at higher temperatures. Additionally, the 
comparatively good agreement between the GSM7 and 
ATOMIC models at the lowest temperature and highest 
density is due to a fortuitous cancellation of the effects as- 
sociated with the approximations employed in the GSM7 
model. 



The behavior observed for the iron spectra is similar 
to that of nickel. While the agreement between the mod- 
els is again excellent up to an electron density of ap- 
proximately 10 19 cm~ 3 (top two rows of Fig. the 
GSM7 results are clearly differing from the other two 
models at densities of 10 21 cm -3 (third row of Fig. [3]) 
for all three temperatures, indicating that the statisti- 
cal treatment has broken down. By an electron density 
of 10 24 cm~ 3 (bottom row of Fig. [3]) all three models 
have diverged for the two lowest temperatures. Much 
like the nickel case, no model containing either the sta- 
tistical treatment or a ground-state-only, quasi-static ap- 
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FIG. 4: (Color online) Ca emission spectra at T e =1.29 x 10 3 , 3.45 x 10 3 , and 6.89 x 10 3 eV (1.5 x 10 7 , 4.0 x 10 7 , and 8.0 x 10 7 K) 
for the three models considered at four different electron densities (top to bottom): 10 10 , 10 17 , 10 19 , and 10 23 cm -3 . The solid 
line (green) was calculated with the ATOMIC model, dashed (blue) with the GSM7 model, and the dot-dot-dashed (red) with 
the GSMF model. 



proximation can be expected to be valid above this den- 
sity for relatively low temperatures. However, as in the 
Ni case, the quasi-static approximation is still valid at 
the highest temperature because the total population of 
the Li-like ionization stage is small and the excited-state 
populations of the He- like and H-likc ionization stages are 
negligible. The comparatively good agreement between 
the GSM7 and ATOMIC models for the lowest temper- 
ature at the highest density is again fortuitous. It is 
worth reiterating that the ground-state-only, quasi-static 
approximation remains valid for a wider density range 
at higher temperatures, underscoring the notion that Li- 



like metastable states are the cause of the breakdown in 
the ground-state-only, quasi-static approximation at the 
lower temperatures. 

While the electron densities for which these behaviors 
are observed in calcium are different, the pattern is once 
again the same. Fig. 0] shows the same sort of excel- 
lent low-density agreement that was described above for 
nickel and iron, but by densities of 10 19 cm -3 (third row) 
the GSM7 model has begun to deviate from the others, 
indicating a breakdown in the statistical methods used. 
By densities of 10 23 cm -3 (bottom row), all three curves 
diverge from each other at the lower temperatures, indi- 
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FIG. 5: (Color online) A comparison of the Ka line ratio 7? predicted by ATOMIC, and the GSM7 and GSMF models as a 
function of electron density for the elements and temperatures considered in this work. The solid line (green) was calculated 
with the ATOMIC model, dashed (blue) with the GSM7 model, and the dot-dot-dashed (red) with the GSMF model. Of 
particular note is the excellent agreement among all three models over the entire range of conditions. 



eating the breakdown of both approximations. Yet again, 
the agreement between the ATOMIC and GSMF mod- 
els at the highest temperature is well understood based 
on the previous explanation, and the comparatively good 
agreement between the GSM7 and ATOMIC models for 
the lowest temperature at this highest density is fortu- 
itous. The pattern that the ground-state-only, quasi- 
static approximation is valid over a greater range in den- 
sity at higher temperatures is also observed in calcium. 

With the comparisons from Figs. (2J-TJ] in hand, it is 
instructive to compare the calculated validity limits with 
estimates provided in previous works. While we are not 
aware of any estimates for the breakdown of the statis- 
tical treatment, the ground-state-only, quasi-static ap- 
proximation has been estimated to be adequate for He- 
like ions up to a critical electron density of approximately 
N e crit = 6 x 10 12 (Z - 2) 4 ' 3 cm" 3 0, where Z is the 
nuclear charge. Beyond this density the Li-like levels aris- 
ing from the ls 2 2p configuration are expected to become 
important contributors via collisional ionization to the 



population of the He-like levels that produce the spec- 
tral lines of interest. This expression gives A*""' 1 * values 
of approximately 7.2 x 10 18 cm~ 3 for Ni, 5.1 x 10 18 cm~ 3 
for Fe, and 1.5 x 10 18 cm -3 for Ca. The models consid- 
ered in this work show that the ground-state-only, quasi- 
static approximation is valid at densities that are one to 
two orders of magnitude higher than those predicted by 
this expression for all three elements at the temperatures 
considered. 

Further inspection of the spectra presented in Figs. 
[2HH indicates that all three models agree in their predic- 
tions of the quenching of the z line as density increases. 
In fact, the values for the density-sensitive line ratio R 
(presented in Fig. O predicted from each of the three 
models are in excellent agreement over a broad range of 
electron densities. The statistical methods used by GSM 
and the ground-state-only, quasi-static approximation do 
not begin to break down until densities for which the R 
ratio has a value of nearly zero. It is worth noting that 
this density range is adequate to model most astrophys- 
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ical and magnetic confinement fusion plasmas. 

V. CONCLUSIONS 

The ground-state-only, quasi-static approximation and 
statistical treatments for excited states arc often-used 
methodologies to simplify the spectral modeling of plas- 
mas at low densities. This work confirms their low- 
density validity, based on detailed atomic models, and 
provides numerical values for the electron densities at 
which these approximations break down for Ni, Fe, and 
Ca at various temperatures near and above the temper- 
ature of maximum abundance for the He-like ionization 
stage. At the temperatures considered, the quasi-static 
approximation using only the ground states of the H-likc 
and Li-like ionization stages appears to be valid for den- 
sities that are significantly higher than those predicted 
by previous works [U, [B[ . The ground-state-only, quasi- 
static approximation is also valid up to densities that are 
orders of magnitude higher than the density limits calcu- 
lated for the statistical treatment employed in the GSM7 
models in all cases considered. This work also confirms 
that the ground-state-only, quasi-static approximation as 
well as the methods employed to treat excited states as 
contributions to effective rates via a statistical treatment 
are valid for predicting the quenching of the z line, and 
thus the drop in the R ratio, as a function of electron 
density in the temperature range considered. However, 
this study further indicates that these approximations 
should be used with care at densities where the R ratio 
would be expected to be close to zero, such as those found 
in laboratory produced sparks (l7j , high-density electro- 
magnetic pinch devices p3l . [25j , and inertial confinement 
fusion experiments (l8j . 

Further study is required to determine if the quasi- 



static approximation using only the ground states of the 
H-like and Li-like ionization stages has a smaller range of 
validity with respect to density than the statistical treat- 
ment at temperatures significantly below the tempera- 
ture of He-like maximum abundance, or for conditions 
found in transient plasmas. We do expect the range of 
validity of the ground-state-only, quasi-static approxima- 
tion to be narrower in density at temperatures lower than 
those considered in this work because the population of 
the Li-like ionization stage should be greater, and thus 
the effect of the Li-like ionization stage on Ka spectral 
formation via collisional ionization would play a more 
significant role. 

Another avenue of future research concerns the im- 
provement of both the statistical treatment and the 
quasi-static approximation used in this work without 
abandoning either approach. The models that employ 
statistical methods could be improved by treating more 
states explicitly and fewer statistically. Similarly, the 
quasi-static approximation employed could be improved 
by allowing for more than just the ground state to be 
populated in the ionization-balance portion of the calcu- 
lation. Investigating both improvements as a function of 
temperature and density is left to future work. 
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